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I. INTRODUCTION 



Physics beyond the standard model (SM) can manifest itself by altering the expected rates 
of Flavor- Changing Neutral-Current (FCNC) interactions. Thus, testing SM and probing 
new physics through top quark FCNC decay is interesting. The top quark FCNC processes 
t — > q + V (y = g, 7, Z) have tiny branching ratios in the SM and are probably unmeasurable 
at the CERN Large Hadron Collider (LHC) and future colliders. Therefore, any positive 
signal of these rare decay events would imply new physics beyond the SM. As the LHC will 
produce abundant top quark events (about 10^ per year), even in the initial low luminosity 
run (~ 10fb~^/year) 8 x 10^ top-quark pairs and 3 x 10^ single top quarks will be produced 
yearly. Thus one may anticipate to discover the first hint of new physics by observing 
anomalous couplings in the top-quark sector. 

Recently, from their measurements of the total cross sections, both DO [l| and CDF ^ 
collaborations at the Fermilab Tevatron have searched for nonstandard-model single top 
quark production and set upper limits on the anomalous FCNC couplings ^f^/A and /tfu/A, 
where the leading order (LO) cross sections have been scaled to next-to-leading order (NLO) 
[sl (or resummation {4]) predictions. At the LHC, ATLAS collaboration has presented its 
sensitivity by studying FCNC top decays These results show that top quark FCNC 
couplings will provide a good probe to new physics beyond the SM. Although there are many 
discussions in the literatures on top quark production and rare decay processes involving 
model- independent top quark FCNC couplings |6l4l8l|. most of them were based on LO 
calculations. However, especially for t q + g, due to the large uncertainties from the 
renormalization scale dependence in its LO prediction through the strong coupling constant 
as, it is necessary to improve the theoretical prediction to NLO in order to match the 
expected experimental accuracy at the LHC. Because of the importance of NLO corrections 
for the experiments I2I, we calculated the NLO QCD corrections to the partial decay widths 
and decay branching ratios of top quark FCNC processes t^q + V{V = g,'-f,Z) more than 
one year ago |l9| . 

In this paper, we describe in detail the calculations in Ref. 19| and consider the FCNC 



operator mixing effects which were ignored in Ref. [19| . This paper is organized as follows. In 
section |Tll we show the relevant dimension-five operators and the corresponding LO results 
in the top sector. Section UlI] is devoted to the detail of the calculation presented in Ref. jiol . 
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Section HVl deals with the evolution of anomalous couplings /tfq. The analytic results for the 
mixing effects can be found in section |Vl while the numerical results are presented in section 
ED 



II. LEADING ORDER RESULTS 

New physics effects involved in top quark FCNC processes can be incorporated in a model- 
independent way, into an effective Lagrangian which includes the dimension-five operators 
as listed below |8[ 

Z 7 

sm Z(7U/ ■^"^ iv ■^"^ i v 

q=u,c q=u,c 

-9s Yl ^^^"'T^ift^ + ^hll^)tGlu + H.c. (1) 

q=u,c 

where A is the new physics scale, 9w is the weak-mixing angle, and are the conventional 
Gell-Mann matrices. The coefficients k^, K^q ^i^d /tfq are normalized to be real and positive, 
while /tq, = 7, (yf) are complex numbers satisfying |/p + = 1. 

From the effective Lagrangian above, we obtain the following LO partial decay width of 
the FCNC top decay in the D = 4 — 2e dimension, 

e 1 



Tlit^q + Z)= .fl [^] (3-/3|-26)-^, 2 
sm 20w \ ^ 1 — e 



where the masses of light quarks q {q = u, c) have been neglected, and Pz = \/l — M'^/m 



a = Here, we also define V^^t ^ q + V) = Vl{t ^ q + \/)Uo, which 



consistent with the results of Refs. 



are 



|8|. 



III. NEXT-TO-LEADING ORDER RESULTS 

Below, we present our calculation in detail for the inclusive decay width of the top quark, 
up to the NLO with the LO partonic process denoted ast q + g. The results of t — g + 7 
and t ^ q + Z are similar to the t ^ q + g, and will be presented all together. 



At the NLO, we need to include both one-loop virtual gluon corrections (Fig. [2]) and 
real gluon emission contribution (Fig. |3]). We use dimensional regularization to regulate 
both ultraviolet (UV) and infrared (IR) (soft and collinear) divergences, with spatial-time 
dimension £) = 4 — 2e. The UV divergences cancel after summing up the contributions 
from the one-loop virtual diagrams and counterterms according to the same convention 
used in Ref. js]. The soft divergences cancel after adding up the virtual and real radiative 
corrections. To cancel collinear singularities, we need to also include contribution induced 
from gluon splitting to light quark pairs at the same order in the QCD coupling. 

The virtual correction oi t ^ q + g contains UV and IR divergences, which has the same 
form as in Ref. jsl and can be written as (the imaginary part is neglected): 

«s ^ 13 17 11 17 2 .r^ A. 

127r V (^iR eiR euv 3 / 

^6zi'^ + ^dzi"^ + ]^5zf + 5Zg^ + 5Z^.j^ Mo, (3) 

where = r(l -|- e)[{47rfi'^)/mfY , and the UV divergences are renormalized by introducing 
counterterms for the wave function of the external fields {bZ\f^ ^ bZ\f^ ^bz'^') and the coupling 
constants [bl 
as in Ref.jSj: 



constants f^Z^s, ^Z^s^/a)- We define these counterterms according to the same convention 



8zf = -^dJ^-^] (—--)-^D. 



27r \ 3 2/ \euv em J Qn euv 
a. f 1 1 



'2 



37r V^uv eiR 
a. „ / 1 2 



bZf = -^DA— + — + 4 
3vr \euv eiR 

«</A = ^r(i + .)(4.r-i-, 

a 1 

5Z,./j, = SZ^zjj, = -^r(l + e)(4vr)^— , (4) 
"1 ovr euv 

where Nf = 5 is the number of the active quark flavors. Here we also present the necessary 
definition 6Zi^ijj^ and SZ^z used in the calculations of the decay modes t — > g + 7 and 
t — )■ q + Z. Combine these together and integrate through phase space, we get the virtual 
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contributions without UV divergences: 

DTT I efn eiR 



-IR 

13 /, 47r/i2 

In — 2 7i? 



-13 In 



mt 



12 In 



2 

mr 



137i 



53 



53 

AT. 



(ln47r - 7£;) 



^2" 



23 



(5) 



All the UV divergences have canceled in T^j^.^., as they must, but the infrared divergent piece 
is still present. We also show the corresponding results for t — g + 7 and t ^ q + Z below, 



virt 



-IR 



eiR 



-2 In 1^ + 27. 



mt 



21ni^ 



In 



mr 



1E 



5 7£; - In 



47r/i^ 



mr 



virt 



37r 



-IR 



eiR 



81n/3z-21n 




27£;-5 



- 21n 



mr 



+8 In 



47r/i^ 



mr 



7s + 



mr 



+4Li2 



/3 



/3| 



-^)ln/3z-ln. ^ 



ll + 5^E-^- 12 
o 



(6) 



The contribution from real gluon emission {t — )■ q+g+g) is denoted as FreaK^ — ^ ?+fi'+fi')- 
In order to cancel all the collinear divergences in the sum of virtual and real radiative 
corrections, we also need to include the contributions from gluon splitting into a pair of 
quark and antiquark in the collinear region, which is denoted as FreaK^ q + q' + q') . Note 
that there are two configurations of final state when the flavor of the light quark coming 
from gluon splitting is the same as the light quark directly from the FCNC coupling, and 
only one configuration when they are different. The contributions of real gluon emission 
(t q + g + g) and gluon splitting {t q + q' + q') are, respectively. 



Tre^x{t^ q + g + g) = —Tl{t-^q + g) 

OTT 



13 

.^IR 



+ - 



13 



47r/i2 
In — 7. 7£ 



mr 



1 

eiR 
53 



13 In 



47r;U^ 



mr 



137i. + f 



In 



47r/i2 



mr 



lE 



31 



171 



-TT 



(7) 



and, 

Frcai(t ^q + q' + q') 



Gtt 



Tl{t^q + 



eiR 



Nf + Nf['^E- In 



47r/i^ 



mr 



1 

12 



(8) 



After adding them together, the total real contributions can be written as 



p9 

real 



TrcaK^ q + g + g) + ^ra^lit ~^ q + q' + Q 



a 



^Tl{t^q + g){-^- — 



+ 



13 



13 
T 



^IR ^iR 

47r/i2 y 53 
In — 7. 7s + — 



-131n— ^ + 137ij 



53' 

Nf 



Nf ( In 



mt 



4:71 fl'^ 

mt 



In — ^ 'jE 



mt 



3l7r2 1025 



(9) 



The corresponding real contributions of t — )■ g + 7 and t — )■ g + Z are also shown below 



real 



^IR ^IR 



^r^(t^g + 7)<| — + — 

, , 47r/i2 
+ I In ^ 



2 In 



47r/i2 



5 I In — ^ 7£; 



mr 



27E + 5 
Ttt^ 52 



real 



3^ 



-81n/3z + 21n^-27s + 5 



mr 



+ 



4(41n/3z + 27s-5) ln/3^ + ln- 



47r/i^ / 47r/i 



mr 



In 



mt 



8 In 27s + 5 +4Li2(/3 



4(l-/3|)(l-4/3| + /3|) 
/3|(3 - 



ln(l - + 



12 + 135/3| - 43/3^ 
3/3|(3-/3|) 



+ 7i? - 57i? 



llTT^ 



(10) 



Combine the real and virtual contributions, we obtain the full NLO corrections for t — ?■ g + (7, 
t — )■ g + 7 and t — j- g + Z as 



rNLo(^ g + 7) 
rNLo(t + 



virt ' real 

a,9 



727r 



1741n( ^ 



wit 



mi 



a 



p7 
real 



-61n( At 



virt 



real 



37r 



r^(t->g + Z) 



-2 In 



mt 



12iV,ln ( ^ 



Att' + 16 



36iV/ - 387r^ + 749 



mt 



4(l-/3|)(l-4/3| + /3|) 



/3i(3 - 



ln(l - 



+4 21n/3^ 



9-/^1 
3-/3I 
12 + 135/31 - 43/3| ^ 2 



ln/3z + 4Li2(/3|)+4Li2 
27r2 - 12 



(3 



/3| 



3/3|(3 - /3|) 

Note that the last expression in Eq. (fTTl) differs from the Eq. (9) in Ref. 19j, and the 
numerical difference is about 2%, as shown below. But the conclusion on branching ratio 



'111 
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for t — )■ g + Z in Ref. 19j is not changed, i.e., the NLO correction is minuscule in branching 
ratio for t — )■ g + Z. Thus, up to the NLO, the partial decay width of the three FCNC decays 
can be obtained 



(12) 



In order to study the effects of NLO corrections to the decay branching ratios, we define the 
following branching ratios for later numerical analysis: 



BR^Lo{t^q + V) 



T{t^q + V) 



(13) 



T{t^W + b)' 

The decay width for the dominant top quark decay mode of t — ?■ + 6 at the tree level and 



the NLO can be found in Ref. 



20|, which we list below for the convenience of the reader. 



8V2' 



T{t^W + b) = To{t ^W + b)<l 



2a. 



/3^)(2/3^-l)(/3^-2) \ 



/3^(3 - 2/32,) 



ln(l - 



9-4/3 



w 



3 - 2/32 



ln/3^ + 2Li2(/3^)-2Li2(l-/3^) 



w 



2/32,(3 - 2/32.) 



TT 



(14) 



where /3vk = ^1 — niy^/ml. 



IV. RENORMALIZATION GROUP EQUATION IMPROVEMENT 



Due to the large scale dependence in the process t ^ q + g, we use the renormalization 
group evolution to improve the result of perturbation theory. The anomalous couplings 
satisfy the following renormalization group equation 

dftS 



dln/i 



(15) 



where 



dg's 



(16) 
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For simplicity, we neglect the subscript tq of the anomalous couplings above and in the 
following discussion. Zi^g^i{gs) is the residue of the renormalization constant 5Z^g. Thus 
according to the Eq. (jlj), we have 

IkO^Qs) = (17) 

Substitute it into Eq. (fT5l) . we can solve the renormalization group equation and get 

= (18) 



where /3o = H — f^/- For cts(/^)) we take it by solving the following renormalization group 



equation 



We do not consider the higher order effects in the 7^9 and /3(as) here because their effects 
are small numerically. 



V. CONTRIBUTIONS TO t ^ g + 7 AND t ^ q + Z FROM FCNC OPERATOR 
MIXING 

The NLO contributions given in Eq. (ITT]) are proportional to the LO results. However, for 
the decay processes t — ?■ g + 7 and t ^ q + Z, there are in general contributions induced from 
the mixing of operators. Because we have no idea about the magnitude of the coefficients 

, these operator mixing contributions may be significant in some cases. In this section we 
will present the contributions coming from the operator mixing so as to complete the full 
0{as) corrections to the decay processes t — )■ g + 7 and t ^ q + Z. 

The contributing Feynman diagrams are shown in Fig. H] and Fig. |5l In previous section, 
we consider only the contributions from Fig. Hl^a) and Fig. [5]^a), (b). In case that all of the 
three are at the same order, the terms proportional to k^k'^ and n^n^ could contribute 
to t — g + 7 and t q + Z with the same significance as the ones we considered before. 
Thus, we will also present these contributions below to investigate the effects from operator 
mixing. 
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For convenience, we introduce the following abbreviations 











^3 

Qf 

Cl 
C2 



2Qfam' 



Ceiftq Aq + ^tq ^tq) ! 




sin 29w cos 6^^^ 

sin 9w 
QfSinOw, 



^(1-/3|)(3 + /3|), 



(20) 



where Q/ is the electric charge quantum number of the fermions, e.g., 2/3 for up- type quark 
and —1/3 for down-type quark. S3 is the third component of the 811(2)^. Notice that we use 
chiral parameters g^iQ^^i instead of fl^ and in the Z-channel for simplicity, which have 
the following relations. 



5'L ~ /tq ^^tq' 



9r = /tq + ^^tq, 



^?lP + |^}j|- = 2. 



i |2 



Denoting F 



virt,mix\ 



7, Z) as contributions to partial decay widths from the sum of 



Fig. 111(b)- (g) and the corresponding counter terms, we obtain 
2a. 



virt,mix 



F^. . 

virt,mix 







37r 



+ 47£; -41n— V - 11 + — 



mr 



25Z, 



97 



4(2 - /3|) 



ro^/3|(3-/3|-26)^^ci 



47r/i^ 



2(l-/3|)ln(l-/3|)-9 




/5|(3 - /3 
+27£;-21n 
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where Cq and the renormahzation constants SZg^ and SZgz are given below: 

_ ^ (1 + /3|)(2 - 2/3| - + ^/3|(2 + ^ + /^l + ^^/^ 

+2Li2(-/3|) + Lis (1 + /3| + + Li2 (^^ (1 + /3| - ^a;/^) 

-Li2 (-1 + 2/3| + /3| + z(l + /3|)a;;3) 

-Li2 (-1 + 2/3| + /3| - ^(1 + 

rv 4 
= -^r(l + e)(47r) 



Svr euv 

fir(l + e)(4vrr^ 
Sir euv 

We can see that T^ij-tmix ^^ai consistent with Eqs. (4.8) and (4.10) in Ref. (2lj, up 



5Zgz = ^r(l + e)(4vr)^Hi±^. (22) 



to a factor of |, which is absorbed into Qf in our formula. Besides, T^iy-^^^iy^ and 5Zgz are 
consistent with Tl-^^^^:^^ and 5Zg^ if we take the limit Mz — > and Qf = —S3 = 1. 

It worths to briefly discuss the renormahzation constant SZgz here. To renormalize the 
operator involving the Z boson, we need to introduce two renormahzation constants: one for 
vector current and another for axial vector current (or left-hand and right-hand). Due to the 
similarity between the structure of the anomalous couplings of 7 and Z bosons in Eq. ([1]), 
we expect that T^^^^^^^^ and T^^^^^^^^ are equal in the limit of Mz — > and Qp = —S3 = 1. 
Thus, we introduce the renormahzation constants at the level of the decay width in contrast 
to the usual practice which is done at the level of the scattering amplitude. Nevertheless, 
this treatment would simplify our calculation, so we take it as the definition of SZgz in this 
paper. 

Denoting J^\{i = 7, Z) as the sum of Fig. [5] (a) and (b), and A^^(« = 7, Z) as the sum 
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of Fig. |5] (c) and (d), we can present the contributions T^.^^^ 
I EM^d$ = rf^{2n^ - 25), 



as 



p7 
real, mix 



real, mix 



+ 21/3| - 72)ci + (24 - 21/3| - /3|)c2] 



Stt I 3 

+2u;^[-(l + 3/3|)ci + (1 - /3|)C2] tan^ 
+ (1 - ln(l - /3|)[(9/3| - 21)ci + (5 - /3|)c2] + 8ci 

X 



tan 



Li2(/3|) + Li2 ( ^ (1 + /3| + ^a;;^)) + Li^ (l + /3| - ^u;/^) 



Combine these results together, we have 



mix 



mix 



■^7 

virt,mix ' ^ real, mix 



+ r 



7 



n^f61n^ + 2vr^-29 



97r 



12- 



r 



virt,mix ' real, mix 



^ro^ { Y H-72 - 42/31 + 11/31) + C2(24 - 48/3| + 11/3| 
+ (1 - /3|)[(2/3| + 9/3| - 21)ci + (5 - /3|)c2] ln(l - /3|) 

^1 [/5|(3 - /3|)C2 - 2(2 - /3|)ci] 



+2/3 



2 '^/S 



■ tan 



1 + 



+2a;4-(l + 3/3|)ci + (l-/3|)c2] 



tan 



tan' 



+2/3|(3 - /3^)(ci + C2) In ^ + 8ciX 



where 



(23) 



(24) 



X = Li2 (-1 + 2/3| + /3| + z(l + /3|)a;^)^ + Li2 (^^ (-1 + 2/3| + /3| - z(l + /3|)a;^) 



-L.2(/3|) - 2L.2(-/3|) - m i±^i±^ m (l + /^l)(2-2/3|-/3|)+./3|(2 + /3|)a;,_ 



(25) 



VI. NUMERICAL ANALYSIS TO TOP FCNC DECAY 



For the numerical calculation of the branching ratios, we take the SM parameters as given 
in Ref. ^ : 

mt = 171.2GeV, Nj = 5, mw = 80.398GeV, 
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TABLE I: Numerical results of branching ratios. Here fi = nit, ^tq/^ — ITeV ^, /tq*/tq + ^7q^tq ~ 
1 and = gi*gi = 1. 

BR BRlo BRnlo BRtot BRnlo/BRlo BRtot/BRLo 

t^q + g 1.0010 1.1964 - 1.195 

t^q + 'j 0.0544 0.0542 0.0486 0.996 0.893 

t-^q + Z 0.04484 0.04480 0.0459 0.999 1.025 



TABLE IL Numerical results of the partial decay width. Here ji = irit, /^J^/A = ITeV ^, ft^ftci + 
= 1 and g|*4 = gfg£ = L 

Width[in unit GeV] Flo T Ftot r/FLo Ttot/rLO 

t^q + g 1.443 1.577 - 1.09 

t^q + 'j 0.078 0.071 0.064 0.91 0.82 

t^q + Z 0.065 0.059 0.061 0.91 0.94 

mz = 9L1876GeV, a = 1/128, sin^ Ow = 0.231, 

Vth = 1, Gf = 1.16637 X 10"^GeV-2. 

We analyze our results by choosing a special set of parameters and fix fi = rrit in the 
following analysis unless specified. Table Hand Table HTl show the 0{as) effects to various 
decay branching ratios and partial decay widths, respectively, where Ftot = T + Fmix and 
BRtot = rtot/r(^ W + b). From Table [H we see that the NLO correction increases the 
LO branching ratio by about 20% for t — )■ q + g, while the NLO corrections are much 
smaller for the other two decay modes. But after including the operator mixing effects, the 
branching ratio for t — > g + 7 can decrease by about 10%, assuming -j^ = = ITeV 
and /tq*/tq + h^*h1^ = 1. For t ^ q + Z , the branching ratio can increase by about 2%, 
assuming -j^ = -j^ = ITeV" and 9^*9% = gL*9L ~ ^- From Table HT] we can see that for 
the partial decay width the NLO results modify the LO results by about 9% in magnitude 
for all the three modes, and the operator mixing effects can decrease the partial width by 
about 9% and increase by about 3% with the above assumptions of parameters for t — )■ g + 7 
and t ^ q + Z, respectively. 

For convenience, we show the branching ratio of t — )■ g+ (7 as a function of k^/A in Fig. El 
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Using the upper limits measured by the DO Collaboration at the Tevatron we get the 
following results: 

< O.lSTeV"^ ^ BR(t c + g) < 2.69 x 10'^, 

g 

^ < 0.037TeV"^ ^ BR{t u + g) < 1.64 x 10^1 (26) 

Using our previous results js], 0, [l9|, CDF Collaboration presented a more precise results for 
the anomalous couplings and the branching ratios in a recent letter 2| 

-j^ < 0.069TeV~^ ^ BR{t ^c + g)<5.7x 10^^ 

< O.OlSTeV"^ BR(t ^ n + ^) < 3.9 x 10"^ (27) 

Following the analysis in Ref. jsl, we plot the coupling nlJK as a function of the branching 
ratio in Fig. [TJ where the ATLAS sensitivities for the two different expected integrated 
luminosities are also exhibited. From Fig. [TJ we can see that the NLO prediction improves 
the sensitivities of the LHC experiments to measuring the top quark FCNC couplings. With 
an integrated luminosity of lOfb"'^, the ATLAS experiment sensitivities can be translated 
to the following relations on FCNC couplings: 

^9 

BR(t ^ g + ^) > 1.3 X 10"^ ^ ^ > 0.033TeV"\ 

7 

BR(t -> g + 7) > 4.1 X 10"^ ^ > 0.028TeV-\ 

BR(t ^ g + Z) > 3.1 X 10"^ ^ ^ > 0.083TeV-\ (28) 

and with an integrated luminosity of lOOfb"^, they can be translated to the more stringent 
relations: 

^9 

BR(t ^ g + ^) > 4.2 X 10"^ ^ ^ > 0.019TeV-\ 
BR(t g + 7) > 1.2 X 10"^ ^ ^ > 0.015TeV-\ 
BR(t ^ g + Z) > 6.1 X 10"^ ^ ^ > 0.037TeV-^ (29) 

Finally, we illustrate the fact that the NLO prediction reduces the theoretical uncertainty 
in its prediction on the decay branching ratios and partial decay widths of the top quark. 
We define i?Lo(/^) = ro(/^)/ro(/x = rrit) and -Rnlo(/^) = r(/i)/r(/i = rrit), and show the value 
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of R{fi) as a function of fi for t q + g in Fig. [HI It shows that the theoretical uncertainty 
from the renormahzation scale dependence can be largely reduced to a couple of percent 
once the NLO calculation is taken into account, so the NLO results give much more reliable 
theoretical predictions. 

To investigate the contributions from the operator mixing effects, we present the contour 
curves for the variables Re{p*f3) (or Re{gl*gL)) and Re{h'^*h^) (or Re{g^*gB)) for the 7 
(or Z) channel, where f = f = 1 TeV^^ (or ^ = f = 1 TeV"^). In Fig. M (a), we 
show the pure operator mixing effects to the branching ratio of t — > g + 7, while in Fig. [9] 
(b) the total result are shown, which includes the LO, NLO and mixing effects all together. 
Similarly, we present the results for t q + Z in Fig. [TUl 

Considering the mixing effects, we need to modify the inequalities involving 7 and Z in 
Eq. ([28]) and Eq. ([29]), which read as 

BR(t ^ g + 7) > 4.1 X 10"^ 

(t") ' - (t") (t") + ^'<^*^'^) > ^ 10-^TeV-^ 

BR{t^q + Z) >3.1 X 10"^^ 

^ j + [-3.6Re(^f ^i) + 6.2Re(^|*^^)] j x 10^^ > 6.9 x IQ-^TeV"^ 

(30) 

for the integrated luminosity of lOfb"^, and 
BR{t ^ g + 7) > 1.2 X 10"^ 



" ' 0-1 I -r I ( 4^ I Mff:ff, + h^^K,) > 2.2 X 10-''TeV-^ 



-5 



BR{t^q + Z) >6.1 X 10 

^ j + [-3.6Re{gt9l) + Q.2Re{gf,*g^^)] (^^ j x 10"^ > 1.4 x lO-^TeV"^ 

(31) 

for the integrated luminosity of lOOfb"^, respectively. 
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VII. CONCLUSIONS 



In order to perform consistent studies for both the top quark production and decay via 
FCNC couphngs, we have calculated the NLO QCD corrections to the three decay modes of 
the top quark induced by model-independent FCNC couplings of dimension-five operators. 
For t —> q + g, the NLO results increase the experimental sensitivity to the anomalous 
couplings. Our results show that the NLO QCD corrections enhance the LO branching 



ratio by about 20%, as presented in Ref. 19|. Moreover, the NLO QCD corrections vastly 



reduce the dependence on the renormalization scale, which leads to increased confidence in 
our theoretical predictions based on these results. For t — )• g + 7 and t — > g + Z, the NLO 
corrections are minuscule in branching ratios, albeit they can decrease the LO widths by 
about 9%. However, if we further consider the effects induced from operator mixing, they 
can either be large or small, and increase or decrease the branching ratios for t — > g + 7 and 
t q + Z, depending on the values of the anomalous couplings {nl^'^ /A, jl^''^ and hf^'^). 
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Note added. After we completed this work we became aware of a paper by Jure Drobnak, 



Svjetlana Fajfer and Jernej F. Kamenik (larXiv: 1004. 0620]) 23| . where they consider the 
same effects from FCNC operator mixing. 
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FIG. 1: Tree-level Feynman diagram for f — >• g + F. Here, we use q to represents the up-quark and 
the charm-quark. 




FIG. 2: One-loop Feynman diagrams for t ^ q + g. Here we use q to represents the up-quark and 
the charm-quark. 

[23] J. Drobnak, S. Fajfer and J. F. Kamenik. la^-Xiv: 1004. 0620) [hep-ph]. 
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FIG. 3: Feynman diagrams of real gluon emission and gluon split. Here we use q to represents the 
up-quark and the charm-quark. 




(a) (6) (c) id) 




FIG. 4: Feynman diagrams of virtual corrections for t — )• g + 7 and t ^ q + Z . Here we use q to 
represent the up-quark and charm-quark. 




(a) ib) (c) (d) 

FIG. 5: Feynman diagrams of real gluon emission for t — >• g + 7 and t ^ q + Z . Here we use q to 
represent the up-quark and charm-quark. 
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FIG. 8: The ratio R as a function of the renormalization scale. Here, ^ = O.OlTeV 




Re(/7*/5) Re(r*/^) 
(a) (b) 

FIG. 9: The contour curves of the decay branching ratio BR(t — )• q + 7) versus the variables 
Re(/'^*/^) and Re(/i'''*/i^). We show in (a) contributions induced from operator mixing effects; (b) 
contributions including LO, NLO and mixing effects all together. Here, we set x ~ X ~ -'■ TeV~^ 
for simplicity. 
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(a) (b) 

FIG. 10: The contour curves of the decay branching ratio BR(t — t- q + Z) versus the variables 
Iie{gf*gf^) and Iie{g^*gpi). We show in (a) contributions induced from operator mixing effects; 
(b) contributions including LO, NLO and mixing effects. Here, we set ^ = x ~ ■'■ TeV~^ for 
simplicity. 
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